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Abstract 

In this paper it is shown that a generahzed circulant matrix under- 
hes every weakly Coupled Map Lattice (CML), independently of the 
form of the coupling term. Therefore, this matrix will appear always 
perturbative methods are used to get the analytical solutions. In fact, 
the inverse of this matrix provides the analytical solution of the CML 
after using first order approximation methods. This inverse matrix, of 
arbitrary order, is explicitly calculated, thus providing the analytical 
expression for the temporal evolution of the CML. 
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1 Introduction 



Since some decades ago chaos represents an authentic revolution in the way 
many problems are confronted. Its tools and applications have spread over 
all sciences from mathematics to physics, passing through chemistry, biol- 
ogy, economy, so on. At the beginning, chaos only studied one-dimensional 
systems. The most remarkable and well-known work of this class is the 
Feigenbaum's job [T], [2] about asymptotic behaviour of unimodal function 
ruling equations like 

Xn+l = f{Xn) (1) 

In spite of the amazing analytical results in one-dimensional systems, 
they are not suitable to describe collective behaviour, because they do not 
take into account interaction among elements of a system, so, they can not 
describe many natural phenomena. One way to solve this problem, and take 
into consideration interaction in a system, is by using Coupled Map Lattices 
(CML) [3]. 

The explicit form of CML can be given by the equations 



X,(n + 1) = (1 - a)f{X,{n)) + - ^ f{X,{n)) z = 1, ■ ■ ■ , m (2) 

^ i=l 

where Xi{n) i = 1, ■ ■ ■ ,m represents the state of elements in the system at 
time n, f rules the individual dynamics of any element Xi{n) according to 
equation ([T]) and the interation among the individual elements is given by 
^ fi^iin)); in other cases a first neighbours interaction is chosen. 
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CML are used to model physical, chemical, economical and biological 
systems [H [5]. From a mathematical point of view, they show many and 
interesting behaviours: travelling waves, synchronized states, period doubling 
cascades, and so on. Nonetheless, most of this behaviours have been observed 
numerically [SI [71 El E] and do not have theoretical support. Furthemore, 
numerical results present another problem: they are restricted by a finite 
and fixed set of parameter values, specific funcions, and a finite number of 
elements in the CML. Therefore, it is unknown whether or not results are 
general results. On the other hand, some models need an arbitrary number 
of elements in the CML to be described properly. 

This problem has been recently solved in [10], where the authors get 
analytical results. To get the analytical solution perturbative methods are 
applied and then a linear system with an arbitrary number of equations has 
to be solved. An arbitrary order functional matrix has to be inverted, at this 
time, the matrix is a generalized circulant matrix, as it would be expected 
and we will show in section "conclusions and discussion". 

The goal of this paper is to calculate the inverse of the generalized circu- 
lant matrix, a matrix that will underlie in every weakly coupled map lattice. 
Futhermore, using the inverse matrix we will explicitly solve the linear sys- 
tem which provides an analytical expression for the time evolution of the 
coupled map lattice. 
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2 Theoretical results 

Given that the CML 

as ^ 

X,{n + 1) = (1 - ea)f{Xi{n)) + — ^ /(^^ ^ = 1 ■ ■ ■ P e « 1 (3) 

P 1=1 

is weakly coupled, we seek a solution in the form: 

i = o,---,p- 1 

where p{x*) = x* i = 1, ■ ■ ■ ,p and P' {x*) ^ 1 (bifurcation point). 
Solving the system ((31) to order e it is obtained: 

a ^ 

-Af'{x*) + Ai+i = ax* + ^i+i ^ = ■ ■ ,P 

P j=i 

whose matricial expression is: 
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The order p matrix of the system 



B 








1 

-fix*,) 1 
-fixl) 1 



V 



















-/'(<) 



(5) 



is not a circulant one, because xl ■ ■ ■ x*. 

The system ^ has p equations and p unknown coefficients whose matrix 
has determinant 

i5i = (-irnm) + (-ir' (6) 



Hence, given that n?=i fi^*) = p' {x*) ^ 1 (by hypothesis), the system 
is compatible and determined for every a 7^ 0. Moreover, as the independent 
term column is not null because x\ ^ x, ^ ■ ■ ■ ^ x^, \i results 



-^2 p Z-i=l -^j 



Xa 



-^4 „ 2^j=\ Xj 








Therefore, the solution of the system (Jll) is different from the trivial one, 
and is obtained by inversion as follows 
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a 



-f'ixl) 1 
-fix*,) 1 








•^2 ^ p 2^j=l Xj 







••• -fix;) ) 



Xl + I Ej=l X* J 

(7) 



The inversion of the matrix (which is not a circulant one) results in the 
following 



A2 

As 



^Ap J 



a 



(_l)P+i(l_(/P(^*)y) 



MN 



(8) 



where the matrix M (as we will show later) is given by 



/ fix*,)... fix;) f'ixl)... fix;) 



/'(xi).../'(x;) 



1 



f'ixl) . . . fix;) f'ixl) fixi) . . . fix;)fixi) 
1 fixi)...fix;)fixi)fix*,) 



V fix*,) . . . fix;_,) fixi) . . . fix;_,) fixi) . . . fix;_,) 
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n<)fixi) 



fixi) fixi)... fix*. 



p-1) J 



(9) 



and by 
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^2 + 1 Ej=l X* 
— rr* _|_ 1 Y^P rr* 



After operating in (I8|) it results in: 



Yj'n=l (( 2;^+n + p 5I]i=l ^i) nf=n /'(^fc+i)) + ( ^fc + p I^f=l ^/*) 

A; = 1, . . . ,p 

(10) 

Not every = 0, because the solution is known to be different from the 
trivial one. 

If some f'{x*) were zero then the determinant of the matrix would be 
(— 1)^+^, according to equation ((6j), and the inverse matrix would still exist. 

We need to show that j^^j^^^M is the inverse of B in order to prove 
that equation (fTOl) is the solution of system ([4]) . Inverting a funcional matrix 
with arbitrary size is quite complicated, it is enough have a look at M (see 
eq. Ml to assume this conclusion. But, if an explicit expression of the inverse 
matrix is known, no matter how it is found, then the easiest way to prove 
that it is really the inverse matrix is by direct multiplication. We use this 
technique in the following theorem. 
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Theorem 1 The matrix D = jj^^jj^jppyyjM , where 



/ fix*)... fix;) f'{xi)...f'{x;) 

1 nx*)...nx;)nxi) f'(xi)...f'(x;)f'(xi) 



M ■- 



f'ixi)...f'(x;)f'(xi)f'(x*) 



V f'{x*)...nx;_,) f'{x*,)...f'{x;_,) f'{xi)...f{x;_^) 



1 

nxDfix'^) 



is the inverse matrix of 



B = 
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-/'(a:^) 1 









1 ••• -fix;) J 

Proof Let 6jj be the elements of the matrix B and m^j the elements of the 

matrix M , then 

-fix*) iii^j 

1 if j = i + 1, iy^p 

1 if z = j = 1 

otherwise 

1 ifj = i-l, z^l 

1 iii=l^ j 

mi_ijf\x*_^) otherwise 
Let us denote by the elements of matrix C, where C = SM, then it 

yields 



m. 



8 



p 

k=l 

hence: 

i) When j = i — it results 

Ci,i_i = hurrii^i^i + 6i,i+im,;+i,i_i = -f'{x*) + /'(x*) = 

ii) When i = 1, j = p it results 

ci,p = 6iimi,p + 6i,2m2,p = -f'{x\) + /'(x*) = 

iii) when z = j it results 

= -/'(^:) nLi + 1 = - nLi + 1 = -in^Dy + 1 

iv) Otherwise 

Cij = burriij + k^i+irrii+ij = -f'{x*)mij + mijf'{x*) = 

Summarizing C = (1 — (/^(x^))')/, with / the order p identity matrix, and 
the theorem is proved. 

Remark. Observe that if, instead of having a global neighbour interaction 
CML, we had any other, the generalised circulant matrix would not change 
because the coupling terms lie in the intependent term column (see equation 
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Conclusions and discussion 



Most of the results in CML are numerical, despite the great efforts to get 
analytical results in this field [IIl[l2j. An interesting theoretical result was 
given by Lemaitre and Chat [13] who proved that global properties emerge 
from local ones in CML. 

This result is very relevant in weakly coupled map lattices because a 
perturbative solution can be seeked from local behaviour. Perturbative tech- 
niques lead to linear problems. In this paper, the theorems are proved that 
allow to get the solution of the CML. 

In the linearized CML given by ^ one would expect to find a circulant 
matrix taking into account the periodic boundary conditions of the CML. 
However, the matrix of system ^ must also take into account the local 
properties, due to the fact that global properties emerge from them. There- 
fore, the matrix of system (111) is one that simultaneously includes periodic 
boundary conditions and local properties, resulting in a generalized circulant 
matrix. As we have pointed out, global properties emerge from local ones, so 
this kind of generalized circulant matrix must underlie every weakly CML. 
Its inverse, explicitly calculated in this paper, gives the sought for analytical 
solutions of the CML, independently of the form of the coupling term (see 
Remark after the proof of the theorem). 
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